Central peak in the pseudogap of high T c superconductors 
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We study the effect of antiferromagnetic (AF) correlations in the three-band Emery model, with 
respect to the experimental situation in weakly underdoped and optimally doped BSCCO. In the 
vicinity of the vH singularity of the conduction band there appears a central peak in the middle 
of a pseudogap, which is in an antiadiabatic regime, insensitive to the time scale of the mechanism 
responsible for the pseudogap. We find a quantum low-temperature regime corresponding to exper- 
iment, in which the pseudogap is created by zero-point motion of the magnons, as opposed to the 
usual semiclassical derivation, where it is due to a divergence of the magnon occupation number. 
Detailed analysis of the spectral functions along the (tt, 0)—(tt, n) line show significant agreement 
with experiment, both qualitative and, in the principal scales, quantitative. The observed slight 
approaching-then-receding of both the wide and narrow peaks with respect to the Fermi energy is 
also reproduced. We conclude that optimally doped BSCCO has a well-developed pseudogap of the 
order of 1000 K. This is only masked by the narrow antiadiabatic peak, which provides a small 
energy scale, unrelated to the AF scale, and primarily controlled by the position of the chemical 
potential. 
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V ■ I. INTRODUCTION 

^ ■ The Fermi surface phenomena in the high-T c cuprates, and especially BSCCO, have been extensively investigated, 
c/2 | and a broad consensus has developed concerning their main features. The Fermi surface is large and hole-like, with 
a simple topology of a rounded square, or barrel, centered around the M-point 1]. Single-particle phenomenology 
is routinely invoked on the ARPES spectra, thus the 'self-energy' and 'damping' are often extracted from the main 
peak as if it were a coherent, weakly perturbed quasiparticle |2j. However, in the underdoped regime and below a 
temperature scale T*, the metallic state is surfeit with low-energy correlations, about whose relevance for cither the 
pseudogap, or the superconducting mechanism itself, there is no general agreement at present. Various experimental 
observations at low energy have been interpreted in terms of stripes 3 , paramagnons Q , phonons [j| , and supercon- 
ducting fluctuations above T c |(|. All these correlations are at present the object of intense scrutiny, mainly with a 
view to ascertaining whether they enhance or suppress superconductivity. 
| Theoretical understanding of the measured electronic spectral functions of high-T c superconductors has received 
J> , significant attention in the context of these efforts 0, S L3- [Hi UM UA UH 13 Oil ■ Physically, conduction occurs in the 
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copper oxide planes, so the most important electronic states are directly accessible to surface probes such as ARPES. 
This naturally allows for a concentration of theoretical effort, especially because the observed spectra offer some 
outstanding puzzles of their own. Such a long-standing issue is the appearance of the pseudogap , observed near 
the vH points for underdoped systems, and its connection with the AF gap at lower doping on the one hand, and with 
the superconducting (SC) gap at lower temperature, on the other. Experimentally, the pseudogap is clearly connected 
with a (n,ir) correlation [15| . and the most natural candidate for its origin are antiferromagnetic fluctuations above 
*^ ■ their transition point Q . 

The present work attempts to connect several aspects of the low-energy phenomenology of the cuprates in the hope 
of realistically constraining the eventual theory of the optimally doped and weakly underdoped state. We adopt an 
effective weak-coupling framework, and concentrate on aspects least sensitive to model details. Our most important 
observation is that the pseudogap does not really disappear at optimal doping, but is instead rather inefficient at 
suppressing part of the spectral strength around the vH points. The unsuppressed strength appears at the Fermi 
level as an 'antiadiabatic' central peak in the middle of a still fairly wide and deep pseudogap. It is the latter 'high- 
energy pseudogap' which indicates the underlying physical scale, while the 'leading edge' scale, associated with the 
central peak, turns out to be incidental to the dynamics. It is primarily controlled by doping. This interpretation 
does not even depend on the pseudogap being due to antifcrromagnctism as such, but only on the fact that the 
dominant perturbing correlation occurs around the wavevector Q = (tt,tt). It is however different than interpreting 
the high-energy 'hump' in terms of bilayer splitting [ill Il7l Il8l | . 
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We do not enter here the important question why the magnetic correlations undergo an essential change at T c . 
Our main aim is to show that when their observed low-temperature behavior is introduced phenomenologically in the 
calculation of the single-electron propagator, the resulting antiadiabatic peak and pseudogap behave consistently with 
the main features of the 'peak-dip- hump' structure, found in experiments on superconducting optimally doped BSCCO. 
In this way our calculation refers to the superconducting state. We only omit the direct effect of superconductivity 
on the single electron propagation, namely the appearance of a superconducting gap. This is justified by the fact that 
the SC gap scale in ARPES is an order of magnitude below the AF scale, manifested by the high-energy 'hump.' In 
order to reproduce typical normal-state ARPES profiles, which do not show a narrow low-energy peak, we only need 
to overdamp the paramagnons. Our work provides a connection between the observed simultaneous appearances of a 
magnetic resonance and of a narrow low-energy peak in the ARPES profile, as the temperature drops below T c . 

Like some other authors [lfl 12(1 l2lj | , here we use an effective weak-coupling (single-band) approach to describe 
the effect of antiferromagnetic correlations on the single-electron propagation. Given that Ud is large in the high-T c 
superconductors, our starting point is the strong coupling limit, and we use the present calculation to develop a 
phcnomenological framework in which the correct physical regime can be identified for the effective weak-coupling 
approach. Section[n]is thus devoted to placing the present work in this wider theoretical context. Section lTTTl describes 
the model results. A comparison with experiment is found in Section llVI Finally, Section [V] is a recapitulation and 
discussion. 



II. THE ELECTRON SELF-ENERGY AND THE CENTRAL PEAK 



A. Separation of charge and spin channels 



We enter a brief theoretical discussion now on the validity of the weak-coupling single-band approach, with a large 
hole-like Fermi surface, when Ud is large. This is the essential input in our calculation, important for its comparison 
with the fc-dependencies measured by ARPES. 

Recently a considerable improvement in understanding the band dispersion measured by ARPES in the high-T c 
superconductors near optimal doping has been achieved by considering the extended Emery model |22j | in the limit of 
large interactions Ud on the Cu-site. The original Emery model of the Cu02 plane is extended by taking into account 
the direct 0-0 hopping t 1 < in addition to the original Cu-0 hopping to and the difference A p d of the O and Cu 
site energies |23[ . 

In the limit of interest \t'\ > t\f A p d this means that the 'broad' oxygen band is weakly hybridized with the Cu level. 
The Emery model then resembles the Anderson lattice model which includes the accurate symmetry of the electron 
(hole) propagation in the Cu02 lattice, either along the 0-0 axis (t') or along the Cu-0 axis (to). Notably, the limit 
\t'\ t^/ A p d, although probably too extreme for physical purposes, corresponds to the Falicov-Kimball model 24], 
also sometimes invoked in the context of high-T c superconductors [23 . l2rl I27L 1281 l2flj . 

The large Ud limit of the Emery model [22] extended by t' was treated in the homogenous mean field approxi- 
mation [23j applied to the slave-boson representation of the Ud = oo Emery model. The usual objection that the 
mean-field slave boson (MFSB) approximation breaks the local gauge invariance required by the slave-boson theory 
was met [23 by emphasizing that the static mean field merely represents the slow component of the slave boson 
field. This latter, allowed by local gauge invariance, only appears as static when particular physical properties are 
calculated, most notably the physical electron band dispersion. Thus the physical dispersion can be represented by 
the usual non-interacting three-band dispersion, but with strongly renormalizcd tight-binding parameters A p f and t 
instead of A p d and to, while t' remains unaffected by the copper on-site repulsion. Most importantly, in this way the 
observed regime A p f w 4|t'| > t naturally replaces the regime A p d > to > t' , inferred from chemical valency analysis 
and high-energy spectroscopy data. The formula for the antibonding electron band is then 

£ (k) = iVpfcos| + V3sin|V (1) 



where 
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P = 12i 2 /i + A8t' 2 f 2 + A 2 



pf : 



Q = lUt'(3t 2 + t'A pf )f 2 - A p/ (18i 2 /i + A 2 



with fi = sin 2 fc^/2 + sin 2 k y /2 and — sin 2 fc a; /2sin 2 k y /2. It is obvious from Eq. (JJJ that the effective near neighbor 
hoppings t and t' enter the dispersion e(k) non- linearly. This is in contrast to those one-band approaches which include 



hoppings to unphysically [31| distant neighbors, but as independent parameters. (Noteworthily, LDA calculations 32] 
show directly that when such long-range hoppings are induced by reduction of multiband to single-band models, the 
effective parameters depend non- linearly on the near- neighbor ones, as in our case.) We use the single band Q from 
the three-band model with this distinction in mind. 

Actually, after taking into account the fast harmonic slave-boson fluctuations around the mean-field saddle point, 
the MFSB band Q decomposes into the narrow resonant band with dispersion e(k) and a spectral density A^, which 
accomodates approximately S holes (doping 5 > 1) on the O-site and one hole localized in the localized state on the 
Cu-site at the energy A p d, deep below e(k) [2i|,|3(j. The observed structure of the resonant band in LSCO is well 
described by the regime \t'\ > tg/A p( j and evolves with doping 6 according to MFSB predictions. 

Once these renormalizations are taken into account, the remaining low-energy analysis concerns only the resonant 
band e(k) containing the Fermi level. In the above slave-boson calculations the spins of the localized holes on the Cu 
sites are taken as paramagnetic. This is justified for large enough dopings 5, when the band-width of the resonant 
band exceeds the magnetic and/or superconducting energy scales, whereas for 5 ~ models of t-J type may be 
more appropriate. Indeed, for optimal dopings the bandwidth is of the order of 1 eV, whereas the magnetic and/or 
superconducting effects occur on scales lower by at least an order of magnitude. This approach is in principle well 
suited to take magnetic energies into account as a perturbation of the main energy scales associated with the resonant 
band. While the formation of the resonant band is associated in the first place with the slow component in the motion 
of the slave boson, its fast component is alone responsible for the weak magnetic couplings. Several calculations of 
this type were carried out before for strongly interacting electron systems. In particular, for the t' = Emery model 
with large Ud, the residual effective couplings were derived explicitly pKjj. The extension of these ideas to finite t', 
in particular to > ip/A^, with a clear distinction between slow and fast components of the slave-boson field is 
currently under way. These residual couplings can be treated, for example, by the perturbational 2D parquet theory. 
At t' = the latter was shown [3J| to lead to ladder-like results ('fast parquet') in most of the space of coupling 
parameters, and to the marginal Fermi liquid only under very special conditions. Here we consider this ladder-like 
regime appropriate for calculating the electron self-energy £ from Fig. ^ where the wavy line represents the spin 
susceptibility \, and the triangular vertex corrections are neglected. 




t I t 

FIG. 1: Phenomenological one-magnon approximation. 

To summarize, for large Ud there is a natural separation of the slave-boson fluctuations into fast and slow compo- 
nents. The latter appear static when calculating the effective dispersion of the electrons, so in fact our mean-field 
slave-boson renormalization of the electronic band parameters corresponds to taking this slow component into account 
in the charge channel [3(| • This explains why the Fermi surface is large and hole-like. The fast component in the 
spin channel is the paramagnon perturbation of the large, hole-like Fermi surface. In the following, we concentrate 
entirely on the latter, neglecting triangular vertex corrections in Fig. ^ 



B. Electron spectral density 

Following the ideas expressed above, the fermion line in Fig. ^ is taken to represent hole propagation in the 
absence of ma gnet ic couplings. In this approach the electron Green's function appearing in Fig. ^ is f ree - It has 
been observed |35l l36| that corrections from self-consistency tend to cancel with vertex corrections in the absence of 
Migdal's theorem, so it appears generally more reliable to take neither into account than only one. Thus the retarded 
Green's function is just 



Gf (k, w) = ^ » ( 2 ) 

R oj - e(k) + fi + nf 



with e(k) from Eq. Q), and the spectral density of the resonant band. Based on the above MFSB considerations, 
we expect a significant k-dependence in this quantity along the Fermi surface, once the strong on-site repulsion is 



explicitly taken into account. In the present work, we concentrate entirely on the vicinity of the vH point, so will 
eventually be absorbed into a coupling constant. 

The wavy line in Fig. ^ is taken to correspond to the simplest form of the magnetic propagator, 



(uj + ij) 2 -uj D (q) 2 

where Q is close to the AF wave vector (tt,tt), 7 is the damping, and up the dispersion 

c^(q) 2 =£5 2 + C 2 |q| 2 . (4) 

Here ui is the band-edge, and c the paramagnon velocity. An upper cutoff loq to the magnons is also introduced, 
corresponding to the extension ujq/c of the magnon anomaly around Q. 

The static magnetic structure factor related to Eqs. (J3J) and (0J is characterized by the value of Q and physical 
correlation length £ = c/uj. This structural factor is measured directly by the elastic neutron scattering [ 3j| |38l or 
indirectly through the nuclear spin relaxation rate Tf" . Both these types of experiment were recently shown [2£| to be 
mutually consistent (in LSCO) when related by the static limit of Eq. (0J. For simplicity however, the incommensurate 
effects (usually associated with 'stripes') will not be included in the present calculation, as they do not seem to be 
important for ARPES results. The value of Q is therefore taken to be (it, ir) in Eq. J3J, and £ is isotropic. 

Turning further to the dynamical features of Eq. J2J), it should be noted that the possibility of a central peak in 
the magnon response, below £3, appearing tog ether with the strong dispersive branches, is not included. We could 
include it fairly easily, following an ansatz 40] slightly different from Eq. J3J). As shown below, such slow correlations 
(usually associated with 'dynamical stripes') are not needed to reproduce the main features of the electronic spectral 
structure observed by ARPES along the (tt, 0)-(tv, tt) line. At fixed c and low temperature, the main parameter tuning 
the AF dynamics in Eq. is the ratio of the damping 7 and the band-edge ui. Magnetic fluctuations are strongly 
overdamped in the normal state, but as soon as superconductivity sets in, a resonance peak appears at 41 meV, 
around optimal doping j^l^. Notably, Morr has obtained the magnetic resonance peak, observed below T c , 
from a mode with the dispersion Q and u> ~ 20 meV, simply by c hanging from overdamping to underdamping. 

The formal expression for the (retarded) fermion self-energy [4^ may then be rewritten as a sum of two terms |55j 

E B (k,w)=-^ ] J du' [ X fl(Q + q^-^)(l-/(^'))IniG^ ) (k-q-Q,c') 

+Gf(k-q-Q,uj-uj')n(uj')ImxR(Q + q,Lo')] , (5) 

where <?k,q is the effective interaction vertex in Fig.^ and Q = (tt, tt). The first term in Eq. J^J) is the boson propagator 
convoluted with the electron response, the second, vice versa. In the high temperature limit the first is negligible, 
because the Bose occupation term dominates the Fermi factor, n(uj') w kT/ui 3> 1 — /(w'). ^ n the low-temperature 
limit kT < ui, which we consider here, both terms may be equally important, with contributions coming from magnon 
zero-point motion. The antiadiabatic central peak is due to the second (boson response) term in both temperature 
limits, as is the lower side wing, corresponding to occupied states. The effect of the first term in the low-temperature 
regime is twofold: it provides the upper wing (unoccupied states) of the pseudogap, and significantly affects the 
position and spectral intensity of the peaks coming from the second term. 

Each of the two terms can itself be expressed as the sum of two contributions, a 'dispersive' part from the propagator 
poles, and a 'diffusive' part from the poles in the occupation factors. The diffusive terms are proportional to the 
damping in the respective response functions. They are not essential for the physics discussed in the present work, 
although we include them in the numerics, when we compare with experiment. The dispersive parts are responsible 
for both the gap and the pseudogap, when it appears. As shown below, the pseudogap Ape can appear not only 
in the high- (kT > uj), but also in the low-temperature (kT < to) regime, while the true gap is always in the high- 
temperature limit, since in the present model it requires ui — > before kT — > 0. The temperature is taken to be 
10 meV, lower than the other parameters in the problem. 

As already emphasized in previous work |40l |43| | . a special physical regime applies in the vicinity of the vH point, 
where the electrons themselves are slow, in fact static at the vH point itself. Then a frequency 'window' appears, 
of the order of the band-edge 2, in which a weakly damped peak survives. This creates an 'antiadiabatic' central 
peak in the middle of the pseudogap, as long as the paramagnon band-edge is finite. For a thorough example of 
the usual adiabatic regime from the side of broken translational symmetry, see Ref. ^4|, while Ref. |4jj describes 
the antiadiabatic regime without translational symmetry breaking. Both of these study charge density waves. In the 
high-T c context, Ref. 0] discusses the high-temperature antiadiabatic case, while Ref. 0] is concerned with shadow- 
band signatures found in the high-temperature overdamped regime (kT >7>(5), which we can also reproduce. The 
latter two references fl9ll43| take the translational symmetry to be unbroken, like the present work. 
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FIG. 2: Thick lines: antiadiabatic central peak in the renormalized electron spectral density A(u>) at the vH singularity 
k = (n,0) with fj, — s v h, for a particularly simple parametrization. Thin lines: k = (n /2, n /2). The circles have the same 
abscissas in all three panels. 

The appearance of the antiadiabatic central peak is first illustrated in the renormalized spectral density A(w) of 
Fig. |21 for a particularly simple parametrization. We put t' = and /i = £ v h, so that the unperturbed system has 
a square nested Fermi surface touching the vH singularity. The product g\ q Ak is taken as constant, independent 
of the position on the Fermi surface. The band-edge u> — 0.03 eV> kT — 0.01 eV~ 7 = 0.015 eV is set to the 
low-temperature regime, as in the rest of the article, while the remaining parameters are irrelevant for the discussion. 
The thin line in the left panel gives the spectral strength at k — (7r/2,7r/2), simply split into two peaks. They have 
quasiparticle signatures, as shown in the right two panels: ReE has a negative slope at the peak positions, and |ImE| 
is small. 

The thick lines show the situation at k = {%, 0), the vH point, for the same parametrization. The maxima in the 
side wings, also denoted by full circles, are evidently incoherent: the corresponding ReE has a positive slope, and 
|ImE| is large. Clearly a central peak has survived at the vH point, protected by the antiadiabatic mechanism. To see 
this, note that when the boson response is peaked around q « in Eq. ©, the main contribution to E at k « (it, 0), 
the vH point, comes from electrons at k — Q w (0, — 7r), the other vH point, where they are slow. 

Thus the central peak consists of vH electrons which do not scatter because they barely move, so for them even 
the slowest available paramagnons are averaged out; this is the antiadiabatic regime. It is clear that it violates the 
Fermi liquid paradigm, since |ImE| ^ at the Fermi level despite ReE = 0. This is because electrons interact with 
dissipative bosons. When the boson damping 7 is zero, |ImE(w)| = for \lo\ < uj at the vH point itself; this was 
checked analytically |43|. It is further important to note that the antiadiabatic peak does not necessarily appear at 
the Fermi level, since it has its own k-dispersion. In the zone, ReE and |ImE| for the antiadiabatic peak behave 
similarly as for a quasiparticle, including the reduced but finite quasiparticle weight. We shall see that the peak can 
reveal its antiadiabatic origin nevertheless, by disappearing with changing k when the electrons involved acquire a 
significant velocity. 

Our calculation can support various physically motivated notions of a pseudogap simply by performing it in different 
regions of parameter space. Because of this, the actual pseudogap obtained in a particular calculation with realistic 
parameters usually behaves as a transitional form between intuitive limiting cases. For example, one might say that 
a pseudogap exists only if the central peak does not cross the Fermi energy. However, that behavior continuously 
transforms into the usual quasiparticle one simply by increasing the energy of the paramagnon band-edge, and in fact 
the central peak can easily cross Ef before the side wings have disappeared. Similarly, the notion that ReE has to 
have a positive slope for a side wing to be called incoherent is less restrictive than requiring |ImE| to be maximal at 
the peak position. Thus the latter notion of incoherence is found as a limiting case, while the former appears across a 
wide range of parametrizations. Finally, one could define the pseudogap by the spectral weight having a minimum at 
the Fermi energy, and a maximum in |ImE| at the same position. The split quasiparticle (thin line) in the left panel 
of Fig. |21 is a limiting case for such a definition, where all the peaks appearing are coherent. However, these peaks do 
not cross the Fermi energy This corresponds to the most common understanding of a pseudogap, where a coherent 
quasiparticle splits in two because of strong scattering at u> = fi, which is the precursor to the new zone boundary 
when uj — * at fixed temperature. 

Let us narrow the usage of the term 'pseudogap' now, to the one found relevant for the present work. We do not call 
the thin line in the left panel of Fig. a pseudogap, in spite of the valley between the two peaks, corresponding to a 
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FIG. 3: Pseudogap with central peak at k = (71-, 0.157r). Thick line: parametrization in the text. Thin line: \i — 0.035 eV. 
Dashed line: kT = 0.1 eV, F = 0.026 eV. Vertical dotted lines, left to right: paramagnon scale (ui — 40 meV), and Fermi 
energy. 

large |ImE|. However, the peaks themselves are 'coherent,' as discussed above. Here we reserve the term 'pseudogap' 
for manifestly incoherent side wings, like the side wings of the thick line, irrespective of the nature or presence of a 
central peak in the middle. Notably, other parametrizations can give three coherent peaks at the vH point, the two 
side ones like at the nodal point, and an antiadiabatic one in the middle; in that case there is no pseudogap at all, in 
the language of this article. 

III. MODEL REGIMES 

At fixed low temperature, the ratio of the paramagnon 'band-edge' ui to the damping 7 becomes the principal 
physical parameter of the self-energy J3J). We shall show below that this number is relevant to account for the 
principal features of the ARPES measurements in BSCCO and YBCO 0] along the (tt, 0)-(-7r, 7t) line, which we 
shall call X-M, in accord with the crystallographic notation for YBCO. We take ui to be 40 meV, in accord with 
experiment, which puts our calculation in the low-temperature limit, kT < lu. This means that the electrons are 
only perturbed by paramagnon zero-point motion. We argue below it is this 'quantum' pseudogap which is actually 
observed in optimally doped BSCCO. The use of the term pseudogap to refer to the destruction of the Fermi liquid 
behavior by mag netic quantum fluctuations is already well established in studies of the antiferromagnetic quantum 
critical point |46j . 

The parameters in expression (J5J are treated semiphenomenologically, i.e. we shall use them primarily to adjust 
the experimentally observed outcomes, but with regard to physically reasonable values. Let us give a standard 
parametrization now, used throughout the article. The renormalized parameters of the Ud — 00 Emery model in the 
hole picture are: copper-oxygen hopping t — 0.3 eV, oxygen-oxygen hopping t' = — 1 eV, and effective copper-oxygen 
energy splitting A p f = 3.6 eV. The Fermi energy is in the electron antibonding band Q of the three-band model [23j. 
The temperature is kT = 0.01 eV, as already mentioned. 

The paramagnon parameters are the band-edge ui = 0.04 eV, damping 7 = 0.015 eV, cutoff uq = 0.15 eV, and 
correlation length £ = c/uj ^3 lattice spacings. The coupling constant is q ^4k = F = 0.077 eV. Its wave-vector 
dependence is neglected, because we concentrate on k « (7r, 0) and q ~ 0. To get a feeling for it, we note that the 
total range of £ around the Fermi surface is ~0.1 eV for this parametrization, which is roughly an order below the 
width of the non-interacting antibonding band QJ. Hence 100 xf/1 eV may conveniently be imagined as 'percent' of 
the non-interacting electronic scale. The chemical potential is /1 = 0.025 eV from the vH point. (Larger /1 means less 
holes.) Individual parameter values are quoted elsewhere in the paper only to denote deviations from the set given 
here. 

The generic form of the pseudogapped spectral function with a central peak is shown in Fig. [3J As long as we 
are in the vicinity of the vH point, some of the spectral strength of the slow electrons survives in the middle of 
the pseudogap, itself of width 2Apq, near the Fermi energy. The persistence of a central peak in the single-loop 
approximation was noted earlier |4tll43l]. There appears an intrinsic 'leading edge' scale, the small distance from the 
central peak to the Fermi energy. When the chemical potential is shifted toward underdoping, this distance increases. 
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FIG. 4: Spectral functions along the X-M line, from the vH point (top curves) to k x = n/a, k y = 0.4n/a (bottom). Left panel: 
uj — 0.06 eV. Middle: high-temperature limit (kT = 0.1 eV, F = 0.026 eV) with Zi = 0.06 eV. Right: same high temperature, 
ui = 0.02 eV. (For high temperature, F must be made smaller, to compensate for the Bose occupation number giving an overall 
factor of kT.) 

The high background observed in ARPES does not appear here. It was obtained in both ARPES 01 an d Raman ^| 
contexts by taking into account the high-frequency slave-boson fluctuations in the charge channel, which we do not 
consider here. (In the latter case (42| this was done in the so-called non-crossing approximation, a somewhat different 
starting point from the mean-field slave-boson one, on which this article is based.) The upper and lower wings at 
Ape > uj may be understood in the semiclassical language as a consequence of the electronic scattering on the nearly 
static, but still not completely ordered AF-like potential induced by the paramagnons. This interpretation implies 
essentially incoherent side wings, with a large Im S and Re £ with a positive slope. (In fact a different structure can 
also appear, with coherent side peaks, as mentioned in the discussion of Fig. El) In the parametrizations used here to 
compare with experiment, the side wings are in fact incoherent, while the physical regime is at low temperature. 

The pseudogap in the high-temperature limit looks quite similar, as shown by the broken line in Fig. EI I n particular 
the energy scale of the side wings is easily adjusted to be the same. The qualitative behavior is however different, and 
since the distinction is important for the phenomenology, we discuss it now. 

In Fig. 01 we show the calculated low-temperature (left panel) and high-temperature (middle and right panels) peaks 
dispersing along the X-M line. We first take the band-edge rather high, uj — 0.06 eV, to emphasize the antiadiabatic 
peak, which picks up most of the spectral strength. Then the only difference between the left and middle panels is 
the temperature, and we note that the antiadiabatic peak is much less dispersive in the low-temperature case, and 
loses strength before crossing the Fermi level. The low-temperature regime is influenced by both terms in Eq. (JSJ) 
equally, with a competition between the Bose and Fermi contributions. In the high-temperature limit the second 
(boson response) term takes over, with a different dynamics, due to the fact that when kT ^> uj > 7, n(uj) in Eq. © 
becomes ss kT/uj, introducing an additional dispersive factor in the denominator. Note, however, that the intensities 
are given without the Fermi occupation factor — we show the renormalized spectral density A(uj), not f(u — h)A(uj). 
This means in particular that the loss of intensity in the left panel is not due to the Fermi surface crossing, in fact 
we shall see (Fig. IHlin the experimental section) that it occurs just as well when the antiadiabatic peak stays away 
from the Fermi surface. 

It is possible to keep the antiadiabatic peak below the Fermi energy in the high-temperature regime as well, by 
lowering the band-edge uj. This is shown in the right panel. Notice that lowering the band-edge in the high temperature 
regime goes toward the opening of a true gap, so the antiadiabatic signal is much smaller, relative to the side wings. 
Otherwise, the dispersion is qualitatively similar to the left panel, especially so when we realize that lowering the 
band-edge flattens the dispersion in the low-temperature kT < uj case as well (as visible in Fig. |SJ), similarly pushing 
the signal away from the Fermi energy. The important qualitative difference in the behavior of spectral strengths is 
however the following: in the left panel of Fig. 0] we notice that the side signal disappears before the antiadiabatic 
one; in the right panel, they disappear together. In fact, the generic behavior in the high-temperature kT > uj limit is 
rather that the antiadiabatic peak disappears sooner. We shall see in the next section that experimental evidence in 
the superconducting state exhibits the low-temperature behavior, providing one piece of evidence that the measured 
response is in the low-temperature quantum regime. 

The second piece of evidence is the effect of the magnon perturbation on the Fermi surface. In the absence of 




FIG. 5: Thin lines: zeroth-order Ud = oo Fermi surfaces. Points connected with thick lines: Fermi surfaces from maxima 
in momentum distribution curves at tv — fi. Dashed lines: zone diagonals. Left: high-temperature pseudogap, kT — 0.1 eV 
and F — 0.026 eV. Right: low-temperature pseudogap. In this figure, fi — 0.015 eV throughout. Points: experimental Fermi 
crossings, large: Ref. |49( . small: Ref. |5Cf . (The value of F in the high-temperature case is adjusted to give a practically 
identical EDC profile at the vH point as in the low-temperature case, in particular the same pseudogap scale.) 



a quasiparticle crossing, the experimental community has developed various alternative criteria to define the Fermi 
surface, one of which is the position of the maximum in the momentum-distribution curve at fixed ui = n |12j . We 
adopt that criterion in Fig. which shows that the effect of magnon perturbation on the zeroth-order Fermi surface is 
qualitatively different for the high- and low-temperature pseudogaps. In the high-temperature (semiclassical) regime 
kT > uj (left panel), the tendency is to change the shape of the Fermi surface so as to follow the zone diagonal in the 
vicinity of the 'hot spots' 0], *- e - the points of intersection of the Fermi surface with the diagonal. Geometrically, this 
means that the angle of the Fermi surface with the zone diagonal decreases. In the low-temperature (quantum) regime 
kT < uj, shown in the right panel, the result is precisely the opposite, the intersection angle increases, and there is 
even a tendency to turn the effective Fermi surface upwards, resulting in a 'flared' shape. While our parametrization 
was chosen for a best fit to energy distribution curves (see below) , and we do not expect detailed agreement with the 
Fermi surface shape, this qualitative difference is of foremost physical importance. It shows, in effect, that 'hot spot' 
scenarios, for example Ref. |5 1| . correspond to the high- temperature regime of Eq. (J3J). They depend on the similarity 
between 'strong' and 'singular' scattering at the nesting wave vector, but this similarity is qualitatively correct only 
when kT > 2, which is not the observed case. The tendency of the Fermi surface to follow the zone diagonal for 
kT > uj is of course a precursor to the diagonal becoming the new zone boundary, when the paramagnons condense. 
As already stressed above, this can only happen in the present model when ui/kT — > 0. While the upturn of the 
Fermi surface has never been clearly observed — there is only one experiment |50j | , the small points Fig. [SJ which 
seems to show such a tendency — the bending to follow the zone diagonal can be excluded with certainty. The Fermi 
surface of optimally doped BSCCO in the vicinity of the vH points is at least a straight line parallel to the T-X line. 
The marked difference between the two final shapes (thick lines) in Fig. |SJ means that in the present model one must 
choose the right-hand panel to reproduce the experimentally observed shape. Thus, both the evolution of ARPES 
spectra in the Brillouin zone and correction to the zeroth order U = oo shape of the Fermi surface in BSCCO seem 
to point in the same direction, that the observed pseudogap is due to zero-point motion of the magnons. 

This discussion can be followed in the form of Eq.|5j In the semiclassical regime kT > uj, the first term is negligible 
with respect to the divergence of the boson occupation number in the second term, which is a precursor to the true 
gap. In the low-temperature regime, the two terms are of the same order. They contribute even at zero temperature, 
because the paramagnon zero-point motion can excite electrons within uj/2 ss 20 meV of the Fermi energy, and since 
the vH singularity is roughly within this range, that means a lot of them. Thus we can violate the conventional Fermi 
liquid picture simply by putting in dissipative paramagnons and letting kT '/uj — > 0, keeping uj at its observed value, 
as discussed on the simple example of Fig. [21 It would of course be interesting to understand why the paramagnon 
resonance should appear when the system goes superconducting. We hope to shed more light on this question when 
we consider the on-site repulsion explicitly, as mentioned in Sect. ILTI 
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FIG. 6: Comparison with experiment on BSCCO (jagged line ^tj), left: T=103 K, right: T=46 K. Calculation (smooth curves), 
left: overdamped magnons (7 = 0.06 eV), right: underdamped magnons (parameters in the text). The scale kT c ~ 100 K is 
approximately given by the double thin cross, which marks the difference in position between the measured and fitted narrow 
peaks. 
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FIG. 7: Experimental log-intensity distribution in optimally doped BSCCO [T^ (left), compared with the present calculation 
(right). The vertical scale is in eV from the Fermi energy, the horizontal one in n/a along the X-M line. White lines: horizontal, 
Fermi energy; curved, unperturbed Ud — 00 dispersion. The intensities are multiplied by a Fermi factor at 100 K. 



IV. COMPARISON WITH EXPERIMENTAL SPECTRA 



In Fig. we show the comparison with experiment, used to establish the parametrization. The jagged lines are 
measured ARPES intensities of optimally doped BSCCO, integrated along the (0, Q)-(ir, 0) line, in the normal 
(T=103 K, left) and superconducting state (T=46 K, right). The smooth lines are both calculated at the (ir, 0) point 
with the parameters in the text, the only difference being the paramagnon damping: left, overdamped (7 = 0.06 eV), 
right, underdamped (7 = 0.015 eV). The difference, marked by crosses, in the positions of the leading peaks in the 
right panel is ^100 K, the superconducting scale. The vertical scales of the two experimental and the two theoretical 
curves are the same, while the relative scale of theory vs. experiment is arbitrary. The single (overdamped) peak in 
the left panel and the side peak in the right panel correspond to a positive slope of Re E and large 1m E, i.e. are both 
incoherent. Only the narrow (leading-edge) peak in the right panel is coherent, with negative slope of ReE and small 
ImE. 

Having fixed all parameters at the (w, 0) point in Fig. [SJ we now compare model predictions of evolution in the 
BZ with a different set of experimental data [l^, nrs t m Fig. The strong non-dispersive structure observed at 
the Fermi level corresponds to the antiadiabatic peak in the calculation. The lower wing of the pseudogap is shifted 
away from the position of the original three-band dispersion, reproducing the experimental high-energy 'hump' scale 
of ^100 meV. The ridge at ~100 meV binding roughly follows the original dispersion, broken and shifted by the 
paramagnon interaction. The antiadiabatic peak is interpolated between it and the upper wing (not visible because 
of the Fermi factor). Notably, the same behavior has been reported along the X-M line in YBCO [isT ]. 

In Fig. |H1 we show the detailed energy distribution curves (EDC's) corresponding to Fig. All the qualitative 
experimental features are correctly reproduced: both the major and minor energy scales, and the downturn (in energy) 
of the antiadiabatic peak as one moves further away from the Fermi crossing. Such an approaching-then-receding of 
the narrow peak with respect to the Fermi level has been noticed in experiment |T^ , and becomes more pronounced 
with underdoping. 

The reference position of the central peak at (n, 0) is at the Fermi energy, as observed in our simple example in 
Fig. El As soon as the chemical potential is shifted, or one looks at other points in the BZ, the peak moves away from 




FIG. 8: Left to right: a) experimental EDC's along the X-M line, from k y = (top) to k y — 0.3657T (bottom), multiplied by 
a Fermi function at 100 K and offset for clarity b) Calculated intensities for the same situation, c) the same for the line 
parallel to X-M at k x = 0.9-7T. d) experiment for the line parallel to X-M at k x — 0.64-7T. EDC's corresponding to the Fermi 
surface crossing from the maximum in the momentum distribution curve are given by a thicker line. The measurement and 
parameters are the same as in Fig. 

the Fermi energy, producing a 'leading-edge' energy scale of the order of the chemical potential, affected of course 
by its own dispersion. It is unrelated either to the primary AF scale, which determines the width of the pseudogap 
by the 'high-energy' side lobes, or the superconducting scale, which does not appear in the calculation at all. Of 
course, as the coupling constant decreases, or band-edge ui increases, the peak begins to turn back into an ordinary 
(weak-coupling) quasiparticle. 

In the right two panels of Fig. [5] we show what happens as one moves towards the T point in cuts parallel to the 
X-M line. We note a significant redistribution of spectral strength, such that the side peak is much stronger at k y = 0, 
the T-X line itself, but quickly loses strength as one moves perpendicularly away from it in the k y direction, parallel 
to the X-M line. Finally at k y = 0.357r/a, only the antiadiabatic peak survives. Experimentally, much the same 
behavior has been observed, with the proviso that it seems to evolve more slowly in the T direction, so the qualitative 
features we calculate around k x — 0.97T here are observed around k x = 0.647T in experiment Again, the wide 

peak has its own approaching-then-receding sequence, similar to the one observed both in optimally doped [p3 | and 
underdoped samples [5^ | . The fact that the calculated qualitative features continue to match closely the experimental 
situation as one moves away from the X-M line into the zone interior, while the quantitative evolution proceeds at a 
different pace, is possibly due to neglecting the k-dependence of the product g\ q Ak in the calculation. 

The principal outcome of the comparison with experiment is that the experimental situation |l2l | in the supercon- 
ducting state corresponds to the model regime of low temperature and low damping kT 7 < uj. This allows us to 
claim that the pseudogapped regime in fact extends to optimal doping. Furthermore, the fact that our underdamped 
curves are calculated in the normal state means that the main qualitative effect of superconductivity on the ARPES 
signal is due to the reduction of magnon damping in the superconducting state. 

V. SUMMARY 

The present work associates the observed spectra around the vH point with the concept of an antiadiabatic central 
peak. It appears in the middle of a pseudogap, representing that part of the spectral strength which is not suppressed 
by the usual adiabatic mechanism of the opening of a pseudogap. The underlying main AF scale Ape is allways that 
of the side wings in Fig. EI which is observed in ARPES as a 'high-energy' feature, or 'hump.' In this way we are 
able to claim that the leading-edge scale, connected with the narrow peak, is not due to any independent physical 
phenomenon. The position of the antiadiabatic peak is sensitive to the chemical potential, which naturally accounts in 
our scheme for the increase of the 'leading edge' pseudogap with underdoping. We can easily recover the conventional 
Fermi liquid by raising the paramagnon band-edge or lowering the coupling constant. At some point the 'hump' scale 
goes to zero, and one recovers the usual weakly perturbed quasiparticle. This is consistent with the observation that 
the Ape ~ T* scale disappears at overdoping, rather than merging with the superconducting scale |53l l54| . 

The pseudogap with an antiadiabatic peak was found in this work to have two physical regimes, low- and high- 
temperature, relative to the paramagnon band-edge. The regime of low temperature corresponds to experiment in 



optimally doped BSCCO. In the model, it gives rise to a dispersion for the antiadiabatic peak which is qualitatively 
different from the bare = oo one, amounting to a non-dispersive 'feature' at a few tens of meV binding energy. 
Thus the observed leading-edge scale need not be entirely due to the superconducting gap. The low-temperature 
regime kT < to is pseudogapped, because a true gap appears when u) — > before kT — > 0, i.e. it is always in the high- 
temperature regime. As long as ui is held fixed, a pseudogap-like situation will occur for kT < u>, without developing 
into a true gap even for the lowest temperatures. However, it may be a true pseudogap, in the sense that the side 
wings are incoherent, and this is the case in the parametrization used here to compare with experiment. Lowering 
the band-edge in the low-temperature regime makes the central peak disappear just like in the high-temperature case, 
thus naturally accounting for the underdoped situation. 

We took much trouble with Figs. 01 and to choose the low-temperature regime, although the magnon mode at 
41 meV is obviously much higher than the temperature. The reason is that our calculation is so simple and generic 
that it fairly represents the perturbation by any bosonic mode which does not have a slow component. There are a 
number of observed low-energy fluctuations, such as stripes, which we do not take into account here. At present, we 
cannot completely exclude a possible role of slow (spin) fluctuations in the electron response measured by ARPES. 
These can be modelled to some extent by a 'central peak' in the boson response, distinct from the dispersive branches 
studied here, and easily introduced in our calculation by a slightly more general parametrization |40| of Im%. However, 
since the observed spectral strength redistribution and Fermi surface shape both correspond to our low-temperature 
regime, we can ascribe the main features of the ARPES response to the dispersive paramagnons. In this way the other 
low-energy phenomena are relegated to a secodary role, possibly having to do with the shape and spectral composition 
of the side wings. Even if the pseudogap were not due to paramagnons, still the conclusion would remain that the 
characteristic energy scale of the relevant boson is higher than the temperature, hence the pseudogap need not of 
itself imply any additional ground-state phenomena. As things stand, we see no reason to depart essentially from the 
natural interpretation in terms of AF paramagnons. 

The neglected quantum fluctuations in the charge channel are also expected to affect the width and shape of the 
high-energy hump observed in ARPES, which appears sharper and narrower in our calculation than in experiment. 
Apart from that, the EDC's obtained along the X-M line have a striking resemblance to experiment in their main 
features. Both the high-energy scale of ^100 meV, and the leading-edge scale of 20-30 meV are correctly reproduced. 
The main intensity pattern, where the peak loses strength further away from the vH point, without ever crossing the 
Fermi energy, is reproduced as well. The intensity shifts between the central peak and side wings are also obtained. 
Finally, the observed slight variation in the narrow peak position, which approaches the Fermi energy and then recedes 
from it, is also found in the calculation. Thus we believe that we have understood the physical origin of the narrow 
lowest-energy signal along the X-M line to be quite general: the electrons giving rise to this signal are slower than 
the perturbing paramagnons, and so escape the adiabatic suppression which opens the pseudogap. We can easily 
recover the actual experimental situation for T> T c , simply by overdamping the paramagnons, which washes out 
the antiadiabatic peak. This scenario is naturally consistent with the fact that a precursor of the narrow peak is 
sometimes observed above T c . Based on the above discussion, we claim that the pseudogap in optimally doped 
BSCCO is in fact fully developed, of the order of 1000 K, and is only masked by the antiadiabatic peak. In this 
way we can view the superconducting correlations as a third scale, an order of magnitude lower than the pseudogap 
'hump' scale, and two orders of magnitude below the on-site repulsion, here taken into account through the overall 
band renormalization. Their interplay with the antiadiabatic leading-edge scale found here, which is of the same order 
of ~10 meV, should be of interest. It remains to be seen whether the physical regime found here to be relevant for 
the low-energy cuprate phenomenology can be consistently obtained from a microscopic strong-coupling approach in 
the presence of an oxygen-oxygen overlap, as outlined in Section ITU 

To conclude, we described the main low-energy signal along the X-M line in optimally doped BSCCO to be due to an 
antiadiabatic central peak. The pseudogap in the same material is an order of magnitude above the superconducting 
scale, and persists below T c . The small separation between the antiadiabatic peak and the Fermi level appears 
naturally in the calculation, even in the absence of explicit superconductivity. It is primarily determined by the value 
of the chemical potential, in a way consistent with its observed variation with doping. 
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